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ABSTRACT 

The Poynting-Robertson Cosmic Battery proposes that the innermost part of the 
accretion disk around a black hole is threaded by a large scale dipolar magnetic field 
generated in situ, and that the return part of the field diffuses outward through the 
accretion disk. This is different from the scenario that the field originates at large 
distances and is carried inward by the accretion flow. In view of the importance of 
large scale magnetic fields in regulating the processes of accretion and outflows, we 
study the stability of the inner edge of a magnetized disk in general relativity when 
the distribution of the magnetic field is the one predicted by the Poynting-Robertson 
Cosmic Battery. We found that as the field grows, the inner edge of the disk gradually 
moves outward. In a fast spinning black hole with a > 0.8M the inner edge moves 
back in towards the black hole horizon as the field grows beyond some threshold value. 
In all cases, the inner part of the disk undergoes a dramatic structural change as the 
field approaches equipartition. 
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1 INTRODUCTION 

Accretion of matter onto astrophysical black holes is a com- 
plex astrophysical process involving several physical sub- 
processes (general relativity, small and large scale magnetic 
fields, turbulence, plasma microphysics, etc.) Accretion pro- 
ceeds in the form of a disk in which material moves slowly in- 
ward until it reaches an inner edge, beyond which it plunges 
in on dynamical timescales. Understanding where the ac- 
cretion disk ends is of paramount importance in the two 
methods for measuring black hole spin, namely fitting the 
thermal X-ray continuum (Zhang, Cui & Chen 1997), and 
fitting the profile of the relativistically broadened iron Ka 
line (Fabian et al 1989; see also McClintock et al. 2011 and 
references therein). 

The inner edge of the disk is expected to be close to 
the radius of the Innermost Stable Circular Orbit (hereafter 
ISCO). It is well known that, for test particles, the position 
of the ISCO depends on the mass and angular momentum of 
the black hole. The situation is more complicated when, in- 
stead of test particles, one considers plasma in orbit around 
the black hole. It has been shown numerically that any tur- 
bulent magnetic field that may thread the plasma modifies 
its dynamics with respect to the idealized case of test parti- 
cles, and the inner edge of the disk moves inward (Penna et 
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al. 2010). It has even been suggested that since the bulk of 
the disk is turbulent, there is no feature in the flow related to 
the ISCO, and the inner edge of the disk moves inward very 
close to the event horizon (Balbus 2012). What is still not 
very clear is what is the role of any large scale dipolar mag- 
netic field present in the immediate vicinity of the central 
black hole. Igumenshchev (2008), Tchekhovskoy, Narayan & 
McKinney (2011) and others have shown numerically that a 
dipolar magnetic field advected inward from large distances 
can grow to equipartition, and can disrupt the process of 
accretion forming a 'magnetically arrested disk' (Narayan, 
Igumenshchev & Abramowicz 2003) interspersed with in- 
ward spiraling bundles of low-density plasma. 

An independent line of thought, the Poynting- 
Robertson Cosmic Battery (hereafter PRCB; Contopou- 
los & Kazanas 1998; Kylafis, Contopoulos, Kazanas & 
Christodoulou 2011), proposes that the central dipolar mag- 
netic field is generated in situ, with the return part of the 
field diffusing outward through the turbulent accretion disk. 
In other words, the origin of the large scale magnetic field 
is situated in the vicinity of the inner edge of the accretion 
disk, and not at large distances as in the latter case. In view 
of the importance of large scale magnetic fields in regulat- 
ing the processes of accretion and outflows, we believe that 
there is still a lot to be learned from the theoretical study of 
axisymmetric accretion within the framework of the PRCB. 
We, therefore, decided to follow on the footsteps of Thorne 
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& Macdonald (1982) and Mobarry & Lovelace (1986) and 
determine the inner edge of a magnetized disk in general 
relativity when the distribution of the large scale magnetic 
field is the one predicted by the PRCB. For simplicity, we 
ignore here the potentially important role of turbulent mag- 
netic stresses (Balbus 2012). We begin in the next section 
with a presentation of our equations and physical simplifi- 
cations of the problem. In § 3 we obtain our first results for 
a non-rotating (Schwarzschild) black hole, and in § 4 we ex- 
tend our study in the most general case of a rotating (Kerr) 
black hole. Our conclusions are drawn in 8 5. 



2 THE RELATIVISTIC MHD EQUATIONS 

Throughout this paper we adopt geometric units where G = 
c = 1. Semicolon stands for covariant derivative, comma for 
partial derivative. Latin indices denote spatial components 
(1 — 3), Greek indices denote space-time components (0 — 3), 
and denotes the spatial part of vectors. In order to derive 
the relativist ic MHD equations for our problem, we follow 
closely the 3+1 formulation of Thorne & Macdonald (1982) 
used by most researchers in the astrophysical community. 
We start with the most general 4-dimensional space-time 
geometry 

ds 2 = g^dx^dx" , (1) 
and the Einstein field equations 

R»u - \g^R = SttT^u (2) 
with the Bianchi identities 

(BT _ Ig^R).^ = , (3) 

where T^ u , the energy- momentum tensor, satisfies 

T ; 7 = . (4) 

We will restrict our analysis to space times where the torsion 
and expansion (but not the shear) vanish. In a cylindrical 
system of coordinates = (t,r,0,<fi) we henceforth also 
assume steady- state, (. . .),t = 0, and axisymmetry, (...),</> — 
0. The general 4-dimensional space-time geometry of Eq. JTJ) 
may be rewritten as 

ds 2 = -a 2 dt 2 + g^idcj) - udt) 2 + g r rd 2 r + geedO 2 . (5) 

Notice that the theory of steady axisymmetric flows in 
ideal General Relativistic Magneto-Hydro-Dynamics (here- 
after GRMHD) around a Schwarzschild black hole in (3+1) 
space-time split formulation was first developed by Mobarry 
& Lovelace (1986). We choose fiducial 'zero-angular mo- 
mentum observers' or ZAMOs who consider hypersurfaces 
t = constant as 'slices of simultaneity'. ZAMOs move with 
4- velocity 

= -(1 ,0 ,0 ,cj) (6) 
a 

orthogonal to hypersurfaces of constant t. They are non- 
inertial observers with acceleration 



U!tU v . 



Here, the lapse function of our observers a = {gt&l ' 94>4> — 
-gt<f>/g<t><t> and gu are in general functions of r 



and 0. Finally, we define the tensor that projects 4- vectors 
onto hypersurfaces of constant t as 

Our hydromagnetic system will be specified by the following 
choice for the energy-momentum tensor 

= pu»u v + ^-{F»F va - -F^F^g^) . (9) 

4:71 4 

Here, = dx^/dr is the fluid 4- velocity as measured by 
ZAMOs, and p is the density of the fluid. Notice that we have 
set here the pressure equal to zero, i.e. we have assumed a 
cold plasma. F^ u is the electromagnetic field tensor which is 
related to the electric and magnetic fields , B^ measured 
by ZAMOs through 

F „u = v » E v _ v v E » + e ^P BxUp m ( 1Q ) 



Here, e^ u \ p = y | det (g^ ) | [fivXp] is the 4-dimensional Levi- 
Civita tensor. The GRMHD equations of motion are ob- 
tained from Eq.j4} supplemented by Maxwell's equations of 
electrodynamics 



F! 



4nr 







(11) 



where J M is the electric current density that satisfies the 
charge conservation condition Jt^ — 0. The electric current 
is related to the electric field through a generalized Ohm's 
law. Eqs. J4) and ([TT]) now read 

V-B = 

-\7 • ( a j) = 
a 

D T B = 

D T e + -\ V • (a 2 S) + W ij *ij = 
or 

D T S + cr • S + ea + - V • (atW) = 
a 

where, 



(12) 



e = T 2 p+^-{E 2 + B 2 ) 

07T 



S = T 2 pv+ -}-E x B 

47T 



W : 



r 2 P v c 



)V+ ■ 



47T 



(E®E + B( 



>B) + ±(E 2 + B 2 )j 



V x (aB) 

47Ttt 



and 



(1 



-1/2 



(1 + 



~2xl/2 
U 1 



9tt) 1/2 , u; 



Here, e is the internal energy, o^ v is the shear tensor of our 
system of reference, v = a~ 1 dx % /dt is the fluid 3- velocity 
as measured by ZAMOs (v = r _1 it), and D T is the Fermi 
derivative (see Appendix for details). We simplify the prob- 
lem even further by assuming a dipolar magnetic field ge- 
ometry, by restricting our analysis to the equatorial plane 
— 7r/2, and by considering only the simplest case with no 
azimuthal magnetic torques. In (r, coordinates, symme- 
try dictates that 

v = (0,0, v*) 

B = (0, B\ 0) with B r i9 / in general , 

E = (£ r ,0,0) with E 6 e ^0 in general. (13) 

For several applications in astrophysics, perfect (infinite) 
conductivity is a valid approximation. In that limit, Ohm's 
law yields the ideal MHD condition 
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= 



(14) 



which is equivalent to the statement that the electric field 
vanishes in a frame comoving with the fluid. 



3 THE SCHWARZSCHILD METRIC 

Before addressing the general case of an accreting spinning 
black hole, we will first consider the Schwarzschild metric 
on the equatorial plane where 



ds 2 



gttdt 2 + g rr dr 2 + 

-a-^V + (i 

r 

+r 2 d6 2 + r 2 d(j) 2 . 



+ 



(15) 



Here, M is the mass of the black hole. The lapse function is 
equal to a = (1 - 2M/r) 1/2 . uo = and = 0. Further- 
more, 



U" = -(1 ,0 ,0 ,0) , 

a 

a 4 = ^(1,0,0) . 

On the equatorial plane, B r — E° 
and according to Eqs. (fTTJl) & ([T4l). 



2 0^ Z7>0 



1 </> o r 
-— V B 

a 



(16) 
(17) 

= due to symmetry, 
(18) 



with B r e ^ and E e ^ 0. The r-component of the last one 
of Eqs. ([T2]) now yields 



1 



2M \ (rv^) 2 M 



r J 



(rB e ) 
4irpr 

(rB 9 ) 2 
4irpr 

= . 



a 2 \r 



a 2 (rv^) 2 



B 



r(rv ), r 



1 + 



(rv^) 



[l - (rv ) 2 



(19) 

Eq. (j!9|) is the generalization of the Newtonian radial force- 
balance in Schwarzschild geometry. Given the distributions 
of B and p with r and in the disk, it yields the distribu- 
tion of the disk azimuthal angular velocity v^. Those distri- 
butions can only be obtained through complex general rela- 
tivistic magneto-hydrodynamic simulations. In this present 
work, however, we are only interested in determining the po- 
sition of the innermost accretion disk layer where the PRCB 
is mostly active. We will thus address this problem analyti- 
cally. 

The vertical thickness h of the inner disk is a major 
uncertain parameter since it depends on several physical 
parameters (thermal pressure, radiation pressure from scat- 
tered photons, magnetic pressure), and is the subject of on- 
going investigation. According to Fig. 1, we can approximate 
the value of the radial component of the magnetic field on 
the upper surface of the inner disk, B r |+^, with the value 
of the vertical component of the magnetic field on the disk 
midplane. Note that B r — on the disk midplane. In other 
words, 



B r 



+h 



-B r 



-rB 9 and B r , 



r 2 B° 



= XrB" 



(20) 




Return field diffusing 
outward through the 
turbulent accretion disk 



/ 



t / >r 

Outgoing large scale dipolar 
magnetic field generated by 
the PRCB in a layer of radial 
thickness 6 




Figure 1. Sketch of innermost disk region where the large scale 
dipolar magnetic field is generated by the PRCB. Solid/dashed 
lines: outgoing/return field respectively. Dotted line: z-axis. The 
PRCB layer is the innermost disk region of radial extent 5. 



We have introduced here the parameter A = r/h. Fur- 
thermore, according to Contopoulos & Kazanas (1998), the 
dipole magnetic field originates in the innermost optically 
thin disk layer of radial extent S such that 



5(jt- 



(21) 



in which photons coming from the central accretion disk 
region penetrate. This is where the PR azimuthal electric 
current is generated. Beyond that layer, the magnetic field 
reverses polarity (see Fig. 1). Here <jt is the Thompson cross 
section of photons scattered by the inner disk plasma orbit- 
ing electrons, and m p is the proton mass. Contopoulos & 
Kazanas (1998) showed that, in order for the PRCB bat- 
tery to operate secularly (i.e. continuously), thus generating 
a large scale dipolar magnetic field interior to the disk in- 
ner edge, the reverse polarity magnetic field needs to diffuse 
outward through the accreting plasma. This occurs natu- 
rally in a turbulent accretion disk with magnetic Prandtl 
number P m ^ 1. We will here assume that 

6~h, (22) 

although this may be generalized. Therefore, in the PR layer, 

-XB° . (23) 



Finally, in order to get rid of the radial derivative of v^, we 
will assume a Keplerian velocity profile, namely 



Under the above assumptions, eq. (|19p becomes 



(24) 



2M\ 



r^2 



+ 



vi(r) 



{(2A-1)( 



-(rv ) 2 (2X 



M 

1 - 
2M 



(A- 



)}[»■ 



(rv^ 
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(25) 



where, v A (r) = (rB e ) 2 /4np. 

In order to determine the position of the disk inner edge, 
we will consider a virtual inward/outward displacement of 
the innermost layer of material of radial and vertical thick- 
ness 5 « h and h respectively, and study its stability. During 
that displacement, mass conservation implies that 



2itr5hp = constant 



(26) 



In order to proceed, we need to make certain further as- 
sumptions about how the various physical quantities in our 
problem vary during the above displacement. Our first as- 
sumption is that the angular momentum parallel to the sym- 
metry axis per unit energy / is conserved, i.e. that 



d0 
dt 



(27) 



Our second assumption is that of flux freezing during the 
displacement. This implies that magnetic flux per unit mass 
is also conserved, i.e. 



= constant. 



ph 



Our final assumption is that 



A = 



constant 



(28) 



(29) 



during the infinitesimal displacement. We see that under 
these assumptions, v\ varies inversely proportionally to r 
during the displacement away from the position of equilib- 
rium, namely 



2 nsco 
v A 



(30) 



where v 2 A = v A (nsco) El 

As we noted before, eq. (|19|) is the generalization of the 
Newtonian radial force-balance, which can also be seen as 
the zeroing of the first radial derivative of an effective po- 
tential V e ff. Therefore, in analogy to the Newtonian case, 
the displacement stability of the innermost disk layer is de- 
termined by the sign of the second radial derivative of V e s. 
Marginal stability corresponds to zero second derivative of 
Veff, i.e. to zero first derivative of eq. ([T9]) . When we substi- 
tute eqs. (|27[) & (|30[) in eq. (|25p . the condition for marginal 
stability, becomes 



T{x, I; A,v A #isco) 
d 



, and 



dx 



T(x, l\ A, va^isco) = , 



(31) 



where, 



T{x, /; A, va^isco 



+ 



Va^isco 

X 2 



{( 2A - 1 )( 1 4) 



)(2A 



2A 



, (32) 



Eq. (I30|) may be generalized as v A (r) « v A (risco/ r ) K 
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Figure 2. Dependence of ccisco = f/M on the inner mag- 
netic field strength for various black hole spin values a. At the 
points the curves end they begin to turn up and backwards. 
Top plot: thin disk (A = r/h = 10). Bottom plot: thick disk 
(A = 1). The various curves from top to bottom correspond to 
a = 0, 0.2, 0.4, 0.6, 0.8, 0.82, 0.84, 0.86, 0.88, 0.9, 0.95, 0.99, 0.999M 
respectively. 



and x = r/M, / = l/M. Notice that in this approach we do 
not need to derive an expression for V r e ff(^)- 

The position xisco of the innermost marginally stable 
circular orbit as a function of our parameters A and v A is 
obtained numerically by eliminating / in the above system of 
equations (see Fig. 2 for a = 0). For < 1, the ISCO grad- 
ually moves outward from its unmagnetized position 6M. In 
other words, the magnetic field acts to destabilize the disk in- 
ner edge, and the disk inner radius must be larger than 6M 
for stability, in accordance with Lovelace et al. (1986). How- 
ever, for higher values of va an ISCO ceases to exist. We 
conclude that, in the presence of a large scale dipolar mag- 
netic field, the ISCO moves outward, and eventually disap- 
pears as the field grows beyond a threshold value. In the case 
of a thin disk (A = 10), that threshold value corresponds 
to va ~ 0.06. We will return to this important point in the 
Discussion section. 



4 THE KERR METRIC 

The above can be directly generalized in the case of a spin- 
ning black hole. In Boyer-Lindquist coordinates, the Kerr 
metric reads: 



ds 2 = gttdt 2 + 2gtd>dtd(j) + gttdr 2 + geedQ 2 + 
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-(1 
E 



2Mr 2 



4Mar 



+ -dr^+Edr + 



E 

^2 



dtd0 



(33) 



where M is again the mass of the black hole, a its angular 
momentum per unit mass (0 ^ a ^ M) and 



; r 2 - 2Mr + a 2 



2.2 2, 

r + a cos / 



A= (r 2 + a 2 ) 2 



2 A 

a A 



The lapse function is equal to a = (AE /A) 
2aMr/A. On the equatorial plane (0 = 7r/2), 

1/2 2 

r 

71' 



1/2 



2M | a 2 



= - ( 1 , , , 

a \ 



2aMr\ 



A 



-) 



Mr 2 



1 + 2- 



,Ma 2 



+ ■ 



,0 ,0 



As before, B r = E = due to symmetry, 



ET = -v*B e 
a 



(34) 
(35) 



(36) 
(37) 
(38) 

(39) 



with B r e ^ and E e e ^ 0. In order to estimate the po- 
sition of the ISCO when the PRCB is active, we will make 
the same simplifying assumptions as in the previous section. 
The r-component of the last one of eqs. (|12|) simplifies then 
considerably and becomes 



-(r^ + wrj 



-%*{£(A+U + A} 



vaHsco 
2r 3 r 2 



{x(v ) 2 + y v + z} = o , 



(40) 



where the detailed expressions for cr r0 , X, y, Z, the 
Christoffel symbols T^, and the flow Lorentz factor V 
can be found in the Appendix. As before, v\ = v^(nsco). 

The assumption that angular momentum parallel to the 
symmetry axis per unit energy / is conserved during a virtual 
infinitesimal displacement of the innermost layer of the disk 
is now generalized as 



= d0 
~~ dt 



-A(a - I) + a(r 2 + a 2 



-a(a - I) A + (r 2 + a 2 ) 2 - la(r 2 + a 2 ) 



(41) 

1 — 

(Bardeen, Press & Teukolsky 1972; hereafter BPT72). When 
we substitute the above in eq. (|40p , we obtain a complicated 
expression of the form 

0(x,Z~;A,vaXisco) = , (42) 
and the condition for marginal stability becomes that 

(43) 



d ~ 2 
— 5(x, /; A, v A xisco) = 



As before, x = r/M, and [ = l/M. The position of the ISCO 
as a function of our two free parameters A and \\ is obtained 



as a solution of the above system of equations (|^2|) & (|l3|) 
for x and /. The solution for various values of the black hole 
spin parameter a is obtained numerically in Fig. 2 for two 
values of A = r/h that correspond to a thin and thick disk 
respectively. Similarly to the case of a Schwarzschild black 
hole, for small values of va, the ISCO always moves outward. 
However, for values of a > 0.8M, an interesting transition 
takes place: as va grows beyond a threshold value, the ISCO 
moves back inward. As we will now see, this result compli- 
cates the observational determination of the black hole spin. 



5 DISCUSSION 

The displacement of the ISCO in the presence of a large scale 
magnetic field generated in situ has direct implications in the 
currently most promising method for measuring black hole 
spin, namely fitting the thermal X-ray continuum (Zhang 
et al. 1997). This method requires knowledge of both the 
Schwarzschild radius and the radius of the inner edge of 
the disk, i.e. it requires knowledge of xisco- As is shown in 
BPT72, in an unmagnetized disk, xisco is one-to-one related 
to the spin of the black hole. In the present work, we showed 
that xisco depends also on the value of the magnetic field 
accumulated in the innermost region of the accretion disk. 
This result introduces one extra complication. 

In particular, the disappearance of the ISCO as the 
field grows towards equipartition (defined in this context 
as va — >• 1) is an indication for a dramatic change in the 
structure of the accretion disk in the presence of a strong 
enough magnetic field generated by the PRCB. It is nat- 
ural to speculate that, when the conditions for marginal 
stability of the inner edge cannot be satisfied, the accumu- 
lated magnetic field will escape outward through magnetic 
Rayleigh- Taylor- type instabilities, as proposed in Igumen- 
shchev (2008), Tchekhovskoy, Narayan & McKinney (2011), 
and Kylafis et al. (2012). As a result, the accumulated mag- 
netic field will decrease, the ISCO will reappear close to 
its unmagnetized disk position, and the disk will undergo a 
new phase of field growth and evolution. For a stellar mass 
black hole, the time required for equipartition field growth 
is on the order of a few hours to a few days (Contopoulos 

6 Kazanas 1998; Kylafis et al. 2012), whereas the intervals 
of dramatic disk structure change are expected to be much 
shorter (dynamic timescales). 

For black holes with spin parameter a < 0.8M, the in- 
ner edge of the disk is expected to be outside its unmagne- 
tized position, therefore the thermal X-ray continuum fitting 
method will on average yield a systematically smaller value 
for the black hole spin. On the other hand, for fast spinning 
black holes with a > 0.8M, the inner edge of the disk moves 
to lower values, and in those cases the thermal X-ray con- 
tinuum fitting method will on average yield systematically 
larger values for the black hole spin. This result introduces 
an artificial segregation in our observational method that 
may be related to the seeming absence of evidence for the 
black hole spin powering the jets in X-ray binaries (Fender, 
Gallo & Russell 2010; see however also Narayan & McClin- 
tock 2012). In other words, the introduction of one extra 
free parameter, namely the amount of magnetic flux ac- 
cumulated around the inner edge of the disk, complicates 
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Figure 3. Range of xisco values as a function of the black hole 
spin parameter a for a thick disk r/h = 1. Dashed curve: unmag- 
netized disk. Solid curve: maximally magnetized inner disk that 
supports an ISCO. Shadded region: range of xisco values when 
^ v^ ^ v A| max - Diffuse cloud: range of ccisco values obtained 
for the X-ray black hole binary LMC X-3 over a period of 25 
years. Left vertical arrow: average spin parameter of unmagne- 
tized disk. Right vertical arrow: spin parameter that yields the 
observed range of xisco values in the context of the PRCB. 

the direct determination of the black hole spin through the 
knowledge of xisco alone. 

The above are summarized in Fig. 3 where we show the 
range of xisco values as a function of a for a thick disk 
(r/h = 1 at its inner edge). The dashed curve corresponds 
to the unmagnetized disk (BPT72), whereas the solid curve 
to a maximally magnetized disk that supports an ISCO. The 
dashed-dotted curve corresponds to the maximum values of 
^isco reached as the inner field (va) grows in a spinning 
black hole with a > 0.8M. In other words, according to 
our simplified model of the inner disk, observations of xisco 
must lie in the shaded region. As an example, in the most 
studied object, LMC X-3, estimates of xisco range from 
about 4 to about 6 (Steiner et al. 2010, assuming a canon- 
ical value for the black hole mass equal to 7.5M©). If the 
PRCB did not operate and the disk were unmagnetized, 
these values would yield an average black hole spin value 
a ~ 0.3M (McClintock et al. 2011). It is natural, though, to 
expect that the PRCB does operate as described above, and 
that the various estimates of xisco simply correspond to 
the various stages of magnetic field growth that the source 
is found to be in. As shown in Fig. 3, in this context, a more 
natural value of the black hole spin in LMC X-3 would then 
be a ~ 0.6M. 

Finally, the reader may wonder in what respect the 
PRCB manifests itself in the inner disks structure differently 
from flux advection from large distances. One way to differ- 
entiate between the two is the radial width 6 of the PRCB 
layer over which B 6 changes significantly. In the PRCB and 
a thin disk, our parameter r/S ~ A is expected to be on 
the order of 10. On the other hand, in the PRCB and a 
thick disk, as well as in a disk where the magnetic field is 
advected inward from large distance, A is expected to be on 



the order of unity. Another way to differentiate between the 
two possibilities has to do with the particular cyclic disk 
variability central to the PRCB, namely intervals of secular 
field growth (hours, days) separated by intervals of strong 
(msec) disk variability. Our present investigation on the role 
of the PRCB in determining the structure of the inner disk 
is certainly not exhaustive and will continue in the frame- 
work of Program ARISTEIA of the General Secretariat for 
Research and Technology of Greece. 
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APPENDIX A: 

Various expressions used throughout the paper are denned 
below according to Thorne & Macdonald (1982): 

V • A = A\ 6 , (V x Ay = e xjk A k , 3 
D T A^ = A^U* - U^U^A^U" 

A - B = A l Bi , (A x BY = e jk AjB k , (A B) ij = A i B j 

Here, e ljk = [ijk]/ ^/|det(7ij)| is the spatial Levi-Civita ten- 
sor. On the equatorial plane [0 = 7r/2), 

A = r 2 - 2Mr + a , S = r 2 , A = r 4 + a 2 r 2 + 2Ma 2 r 
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